Parity-relevant Zitterbewegung and quantum simulation by a single trapped ion 
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Zitterbewegung (ZB), the trembling of free relativistic electrons in a vacuum could be simulated 
by a single trapped ion. We focus on the variations of ZB under different parity conditions and find 
no ZB in the case of odd or even parity. ZB occurs only for admixture of the odd and even parity 
states. We also show the similar role played by the parity operator for the trapped ion in Fock-state 
representation and the space inversion operator for a realistic relativistic electron. Although the ZB 
effect is invisible in a relativistic electron, preparation of the trapped ion in different parity states is 
a sophisticated job, which makes it possible to observe the parity relevant ZB effects with currently 
available techniques. 
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Since the discovery by Schrodinger [T[ , Zitterbewegung 
(ZB), i.e., the trembling of free relativistic electrons, has 
drawn more and more attention and interests over past 
years [2 - 8] . It has been generally believed that the trem- 
bling of a relativistic electron is resulted by the interfer- 
ence between negative and positive energy components, 
a typical relativistic feature of the Dirac electron. Up 
to now, however, no direct observation of ZB has been 
achieved due to inaccessibility with current experimen- 
tal techniques, which led to some questioning on the ZB 
. Moreover, there have been alternative explanations 
for the origin of the ZB, such as the continuously vir- 
tual transition process between different internal states 
in view of quantum field theory Q or the relevance to the 
complex phase factor in context of space-time algebra Q . 

By quantum simulation, some relativistic effects have 
recently been demonstrated in some controllable physi- 
cal systems, such as graphene, semiconductor, supercon- 
ductor and trapped ion |llj - tl8l |. It was recently consid- 
ered that the ZB occurs not only under the relativistic 
condition, but also extensively in the dynamics of a sys- 
tem with more than one degree of freedom [l4|, [HI ■ The 
demonstration of the ZB beyond the relativistic electron 
helps us to further understand the Dirac equation and 
the relevant relativistic phenomena. 

In this work, we focus on the role of parity played in the 
simulation of Dirac equation for the ZB effect by a single 
trapped ion. Like in [TollT6| . we also employ the motional 
degrees of freedom of the ion to simulate the position and 
momentum of the relativistic electron, and the internal 
degrees of freedom to refer to the energy states. Since 
the motional state of the ion could be quantized, we may 
discuss the problem in number-state representation. The 
key point of our work is to introduce a parity operator fl, 
by which we show that, besides the conventional consid- 



eration of the origin of the ZB, i.e., the interference be- 
tween the positive- and negative-energy components, the 
ZB is also relevant to parity of the states. To understand 
the relevant physics, we will compare the parity operator 
of the trapped ion with the space inversion operator of 
the realistic relativistic electron. The experimental fea- 
sibility to observe the parity-relevant ZB effects will be 
justified. 

Under the radiation of three laser lights with red- 
detuning, blue-detuning, and carrier transition, respec- 
tively, the interacting Hamiltonian of a single trapped 
ion reads HH 



H = ihrjVl(a + — a)<r x + Kla z 



(1) 



where r\ is the Lamb-Dicke parameter, f2 and f2 are, re- 
spectively, the effective Rabi frequency and the effec- 
tive Larmor frequency of the ion. a + (a) is the cre- 
ation (annihilation) operator of the quantized motion of 
the ion. a x and a z are usual Pauli operators. Defining 
p = ih(a + — a)/2A with A the size of the ground state 
wave function, we may rewrite the Hamiltonian as a form 
analogous to a f + I dimensional Dirac equation 



Hd = 2ijAQpa x + Mla z 

provided c := 2ifAfl and mc 2 := hfl. 
Introducing the parity operator 



n 



(2) 



(3) 
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commuting with the Hamiltonian in Eq. (JTJ, which 
means n is a conserved quantity under the Hamiltonian 
H, we study below the dynamics of the trapped ion with 
different parity states, i.e., of the odd or even parity, or 
of admixture of the both. 

The ZB is related to the average position (x(t)) of the 
trapped ion. Since our interest is in the states under some 
parity conditions, we have to find the common eigenfunc- 
tions of H and II. To this end, we define \p) as the eigen- 
function of (a + — a) by assuming (a + — a) \p) = —ip \p). 
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Straightforward deduction yields 
1 



\p) 



/2n 



E — 7=r- 1") . 



71=0 



(4) 



where H n (x) is a Hcrmitc polynomial defined by H n (x) = 

^l)n e x 2 /2_£_ e - X 2 /2 ^ j t CQuld bc proven that 

the states in {\p}} arc orthonormal and complete, and 
each state \p) corresponds to the momentum p mom = 
hp I (2 A). 

Substituting Eq.fP into Eq.([T) and diagonalizing H 
yields the eigenvalues 



e± = ±\J{hVL) 2 + {rjhnp) 2 

and the eigenstates 

TjKlp 



\S E+ {p))=N{p) 
\S E _(p)) = N(p) 



1. 



e+\ + nn 
-j]hfip 



\E-\ + m' 



with N(p) = \/(\E± \ + Ml)/ 2 \E±\ the normalization 
factor. So the total eigenfunctions of H are \^e + {p)) = 
(p)) ® |p) and = |5 B _(p)> ® |P>- 

Under the parity operator II, we assume following par- 
ity states in the number-state representation, 

\o) = |+) /a™+i l 2m + !> + I") E / 2 ™ l 2m ) > ( 5 ) 

i e > = i+) E > 2m ) + 1-) E ' 2m + x > ' ( 6 ) 

m m 

where /& (fc = 2m, 2?n + 1, m = 0,1,2,- ••) are nor- 
malized coefficients. |o) stands for odd parity state 
due to II|o) = — |o) and |e) for even parity state with 
II |e) = |e). |±) are the eigenstates of u z with eigenval- 
ues ±1, respectively. 

For our purpose, we may construct the spin states |±) 
by positive energy eigenstates as 

(7) 
(8) 



|+> = W(\S E+ (p)) + \S E+ (-p))), 
\-)=W'QS E+ (p))-\S E+ (-p))), 
or by negative energy eigenstates as, 

\+)=W'(\S E -(-p))-\Se-(p))), 
\-) = W(\S E -(p)) + \Se-(-p))), 



(9) 
(10) 



with normalization factors W = y/\E±\ /(2 \E±\ + 2MI) 
and W -- 



[\E±\ (\E±\ +Ml)]/2(riMlP) 2 . Moreover, 
the odd and even motional states are associated with the 
momentum eigenstates \p), i.e., 

J2f2m+i\2m+l)K-^=(\p)-\-p)), (11) 

m * 

^/ 2m |2m)oc^(b) + |-p)). (12) 



Substituting Eqs. (7), (JSJ, ([TTJ» and (IS) into Eq. ©, we 
may write down the co-eigenstate of E + and odd parity 
to be, 

V4+) = w( P )(\s E+ (p)) + \S e+ (-p))) ® ^(|p) - |-p» 

+QW'(p)(|&») - |& + (-p)» (8 (^(b) + |-p» 



where Q is a coefficient to be determined. We only keep 
some reasonable terms by setting Q = W(p)/W'(p), 
i.e., elimination of the terms of \Se + {p)) \—p) and 
\Se + (— p)) |p)- Then we have, 



Ve+ 



* [|^ + (p))_|^. + (_p))]. (13) 



Similarly, we have other co-eigenstates of E + and even 
parity, E_ and different parities [2lj |. 
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|^| + ) = [|Vb+(p)> + I W~p)>] 



(14) 



= ^= [|^-(p)> + \1>b-(-p))] , (15) 

= ^ [|Vb-(p)> - lV>£-(-p)>] • (16) 

Based on these eigenstates, the average position of 
the ion (x(t)) = ((a+ + a) A) can be calculated by the 
evolved states with odd or even parity, where 



e R 



e R 



|^(i)) = ]Ta p k + (p))e"^ +^6 P U|_(p)) 



with a p and 6 P the coefficients determined by the initial 
condition. 

To calculate (x(t)), we first consider (dx/dt), the aver- 
age velocity, 



dx\ _ A 
dt / ih 



a + + d.H 



D 



2??AO (a x ) 



Taking the odd parity as an example, we have 



E °p 6 p 



E & p a P e " 



2 'i E ±i« 



^_(t) a, r E+ (t)).(17) 
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It is easy to check in Eq. (17) that (i/j°(t)\ d x \ip°(t)) = 
since every term is zero. As a result, (x(t)} remains un- 
changed. Similarly we have (ip e (t) \ d x \ip e (t)) = 0. These 
results imply that a trapped ion in an eigenstate of II 
would be on average static. 

Because Eq. (fl~7|) involves both the positive and nega- 
tive energy components, the ZB should occur, according 
to the conventional viewpoint, due to their interference. 
Our result, however, presents that the ZB depends not 
only on the interference between the positive and nega- 
tive energy components, but also on parity of the states. 

To be more clarified, we have numerically calculated in 
Fig. [1] the average position of the trapped ion with the 
initial state (cos/3 |+) + sin/3 |— )) |0). Since both |+) \n) 
and |— ) |n) are states with definite parity, the case with 
a = /3 is of the most mixed parity. By changing the 
values of /3, we may see clearly from the figure that the 
ZB occurs only in the admixture of the odd and even 
parity eigenstates. 

Our result could be understood by the viewpoint in (l9j 
that the ZB could appear in any system, besides the rel- 
ativistic system, with more than one degree of freedom. 
In our case, the ZB effect is originated from the interplay 
between the internal and motional degrees of freedom 
of the ion. With respect to the conventional viewpoint 
of interference between the positive and negative energy 
components, we may check Eq. (17) again which involves 
both positive and negative energy states. The internal- 
motional-state interplay leads to the interference between 
different energy component terms. Once the ion is in a 
certain parity state, however, the interference is destruc- 
tive, yielding (x) = 0. So the ZB appears only in the case 
of admixture of different parity states, which allows the 
constructive interference between different energy com- 
ponents. Simply speaking, whether the ZB appears or 
not, is decided by both the interference and symmetry, 
the latter of which is reflected by parity. 

For a realistic relativistic electron, there is no demand 
to quantize the motional freedom, but the parity opera- 
tor II discussed above reminds us of the space inversion 
operator P defined as [22j . 



P<p(x,t) = & z lfi(~X,t), 



(18) 



where tp(x,t) is the wavefunction of the relativistic elec- 
tron, and we have used the overline to represent the pa- 
rameters of the relative electron in order to distinguish 
from the ones of the trapped ion. For a relativistic elec- 
tron, the eigenstates of Hd are, 



<PE + (p\x,t) = n (p) 



ip E _ (p;x,t) = N(jp) 



1 

cp 



I E+ 1 +mc' 

— cp 
E-\+mc 2 
1 



e T 



e R 



where m is the mass of the electron, p and x are, respec- 
tively, the momentum and the position, E± stand for the 
energies and N(p) the normalization factor. 
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FIG. 1: (Color online) Numerical simulation of (x(t)) in units 
of A for initial state (cos /3|+) + sin /3| — ))|0), where |±) are 
eigenstates of a~ and |0) is the ground motional state. A c = 
277f2A/$! is the Compton wavelength. We consider (a) A c = 
0.6A, near the non-relativistic limit and (b) A c = 5.4A close 
to the relativistic limit. From the bottom to top, the curves 
correspond to /3 = or tt/2, 7r/12, 7t/6 and 7r/4, respectively. 



We have P(pg (p;x,t) = ±ipjj (—p;x,t). Under such 
a consideration, we may easily find the odd parity states 

<p°e_ = ~^^e- (P) + <pe- i-P)], 

and the even parity states 

V| + = -^PE+iP) + <PE+(-P% 
<Pe_ = -j^VE- (P) - Ve- (-£)]• 

They are formally consistent with the co-eigenstates Eqs. 
(|T3| - (fT6)) for the trapped ion. In this sense, we consider 
that n and P play the same role in the respective sys- 
tems. In other words, n is something like an inversion 
parity operator in the number state representation, which 
moves population back and forth between the internal 
and motional states. 

Following the same steps as for the trapped ion, we can 
immediately find (x) = of the relativistic electron under 
definite parity states. So the ZB of the relativistic elec- 
tron occurs when the electron is in both the admixture 
of different parity states and the admixture of different 
energy components. A particle, such as a relativistic elec- 
tron or a trapped ion, staying in a certain parity state will 
be static or in harmonic oscillation. Otherwise, a moving 
particle with different energy components coexisting will 
surely experience the ZB. 

Since the ZB in a real relativistic electron is inacces- 
sible with current technique, we have to resort to other 
systems, such as the trapped ion, to observe the parity- 
relevant ZB effect. With different trapped ion's motional 
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states, the eigenstates of the parity operator have abun- 
dant forms, such as |±) <g) |n) and |+) (8 \n) A ± h^) £g> \n) B 
where \n) A ( B j is the displaced coherent state [2_|. With 
currently available ion-trap technique, it has already been 
achieved the cooling of the single ion to the motional 
ground state (2^. Preparation of the ion to a certain 
motional Fock state or a certain coherent state has also 
been a sophisticated job [25[. As a result, using the par- 
ity states mentioned above, we may check the variation 
of the ZB with respect to different parity conditions fol- 
lowing the operations in (l6| . 

Besides the simulation of Dirac equation, the parity 
operator n also has application in quantum computing. 
II commutes not only with H, but also with other laser- 
ion interaction Hamiltonians, such as H r = er+a + er_a + , 
a usually used Hamiltonian for logic gate operation (26| . 
If we prepare the ion in the state |+) |0) or other parity 
states, the ion will be more steady during the operation 



than in any mixed parity states. This is advantageous to 
quantum gate operation. 

In summary, we have investigated the dynamics of a 
single trapped ion under some certain parity conditions. 
Our study has shown that the ZB in the trapped ion is 
relevant not only to the interference between different en- 
ergy components, but also to parity. To understand the 
physics related to the realistic Dirac particle, we have 
discussed the correspondence of the parity between a rel- 
ativistic electron and a trapped ion. Experimental feasi- 
bility of observing our results by a trapped ion has been 
justified with currently available techniques. We argue 
that our study is not only helpful to explore the ZB effect 
itself, but also useful to further understand the quantum 
characteristic of the ultracold trapped ion. 
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